


^**tit^ *ii&UD 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

INVESTIGATION OF THE USE 
OF LIQUID CRYSTAL THERMOGRAPHY 
TO STUDY FLOW OVER TURBOMACHINERY BLADES 

by 

Roy L . Brennon 
September 1980 

Thesis Advisor: Matthew D. Kelleher 



Approved for public release; distribution unlimited 



T197818 



i 

n 





UNCLASSIFIED 

seCUMITV CU ASSiriCATIOM or THIS Bum tnfw^) 



NAVJI h: i n ycHui... 

WOwrtRtY CAli^ ^3940 



REPORT DOCUMENTATION PAGE 


READ instructions 
BEFORE COMPLETING FORM 


. NCPOWT mumSCA 2 . OOVT ACCCtSION MO, 


1. AtClAitNT’l catalog numbea 


4. Tl TL C 

Investigation of the Use of Liquid CrysteQ 
Thermography to Study Flow Over Turbo- 
machinery Blades 


*• TYAE op AEAOAT 4 AEAIOO COVEAEO 

Master's Thesis; 
September 1980 


• . ACAPOAMINO OAG. AEAOAT mumBCA 


7 . AuTHOAf»j 

Roy L .. Erennon 


i. CONTAACT OA GAANT MUMBEAr«i 


S. ACAPOAMIMO OAQAMI 2 ATIOM NAME AnO aOOACSS 

Naval Postgraduate School 
Monterey, California 9394Q 


10. AAOGAAM clement. AAOJECT. T ASK 

AAEA 4 «OAK UNIT NUMBEAS 


11 COM TftOI-LIMG O^Yice NAME AMO AOOMCSS 

Naval Postgraduate School 
Monterey, California 93940 


IJ. .e^O.T GATE 

September 1980 


IS. NUMBEA OP A AGES 

67 pages 


fi. MOMITOAimG aGCnCV mmmC 4 AOQACSS/’lf tmm Conlfoillng Office) 

Naval Postgraduate School 
Monterey, California 93940 


1i. SECUAITY CLASS, (oi r»>erfj 

Unclassified 


iS«. OECLAltiPiCATlON/OOANGAAOlMG 

schedule 


l«. OISTHISUTIOM STATCMCmT (oI Ihl* 

Approved for public release; distribution unlimited 


17. OISTAiAuTiOn STATCWCmT fot Bi^ek 30, U kUUrmmt tfmm H 0 p 0 tt) 


It. SUAALCMCmTAAY motes 


IS. KEY WQAOt ('Cdwrinu# «•* fwmf 9*km It nmc999mty kf ktoek n%Mmkkt) 

liquid crystals, thermography, isotherms, stall angle, separation, 
laminar separation bubble, convection heat transfer coefficient, 
boundary layer, angle of incidence, NACA series 65 compressor blade. 


30. abITAaCT (Contlmt^ «*i fmvwrmm mikm It n«c««4AiT AAil ky kikek mmtkkf) 

The use of liquid crystal thermography was investigated as a 
technique for visualizing the flow over a NACA series 65 compressor 
blade. The demonstration of the feasibility of the technique was 
conducted in the low turbulence wind tunnel at the U.S. Naval 
Postgraduate School. 

Local heat transfer coefficients were obtained for Reynolds 
numbers varying from 100,000 to 600,000 with the angle of 



DO 1473 coition OC I MOV «l IS OaiOLCTC 

(Page 1) 0 103-0 14- AAO I : 



UNCLASSIFIED 

ttCU*ITY CLAttIPICATION OA TmU ^AOC ('»»<•" Dtf 



UNCLASSIFIED 




9 



mmm 









Approved for public release; distribution unlimited 



Investigation of the Use of Liquid Crystal Thermography 
Study Flow Over Turbomachinery Blades 



by 



Roy L. Brennon 

Lieutenant Commander*) United States Navy 
B.S. , United States Naval Academy, 1968 



Submitted in partial fulfillment of the 
requirement for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
September 1980 



f w 




ABSTRACT 



The use of liquid crystal thermography was investigated 
as a technique for visualizing the flow over a NACA series 
65 compressor blade. The demonstration of the feasibility 
of the technique was conducted in the low turbulence wind 
tunnel at the U.S. Naval Postgraduate School. 

Local heat transfer coefficients were obtained for 
Reynolds numbers varying from 100,000 to 600,000 with the 
angle of incidence of the blade varying from 0 degrees to 
30 degrees. 
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I. INTRODUCTION 



The two essential turbomachinery parts of the gas 
turbine are the compressor and the turbine. A large part 
of the work produced by the turbine unit goes into driving 
of the compressor. The efficiency of the gas turbine 
engine will decrease dramatically with ah increase in the 
aerodynamic losses in the compressor. 

Aerodynamic losses are irreversibilities in the flow 
caused by profile friction losses, flow separation losses, 
stall losses, and other secondary losses. Therefore, it 
is absoluately necessary to have a clear and thorough 
understanding of the aerodynamic phenomena occurring in 
the compressor blade system. 

The understanding of the flow field through the com- 
pressor allows the blade designer to better compute the 
forces acting on the blades, the location of potential 
aerodynamic problem areas, and determine favorable blade 
shapes with lower losses. Chapters 16 and 17 of Ref. 1, 
contain a complete discussion on the necessity of the de- 
signer to understand the flow field over a blade, as part 
of the design of turbine and compressor units. 

The puirpose of this investigation was to demonstrate 
the feasibility of using liquid crystal thermography to 
visualize the surface characteristics of the flow over 
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a NACA series 65 compressor blade placed in a forced con- 
vection environment. 

Cholesteric liquid crystals, when applied to a black 
surf ace, exhibit brilliant changes in color over discrete, 
reproducible temperature bands. Used as a temperature 
sensor, this allows one to visually observe selected iso- 
therms on the surface. These isotherms can be used to 
infer the location of points of flow separation and 
boundary layer reattachment. The area of flow transition 
and angle of stall could also be determined. 

As documented in Ref . 2 , earlier experiments using 
liquid crystals in wind tunnels to study heat transfer 
rates showed promise but were limited in many ways. The 
pure liquid crystals deteriorated rapidly with age and the 
accuracy was strongly influenced by mechanical shear, 
viewing angle, ultraviolet light, and chemical contamina- 
tion. An encapsulating process developed by the National 
Cash Register Company (NCR) using a gelatin in a polyvinyl 
alcohol binder coating greatly reduced or eliminated these 
problems. These crystals, referred to as encapsulated 
liquid crystals, are small spheroids with diameters on the 
order of 20-50 microns. 

Previous work at the Naval Postgraduate School using 
encapsulated liquid crystals included an investigation of 
the temperature field produced by a cryosurgical probe 
embedded in a clear gelatin-water test medium [Ref. 3] . 
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Experiments using heated cylinders, in a crossflow of air, 
covered with the encapsulated liquid crystals were mapped 
by Meyer in Ref. 2 and Field in Ref. 5. A thesis 
by Durao [Ref. 6] , investigated heat transfer in straight 
and curved rectangular ducts using liquid crystal thermography. 

In the present investigation, data were obtained on a 
NACA series 65 compressor blade covered with carbon impregna- 
ted paper that had a uniform electrical resistivity. The 
surface of the blade was electrically heated by passing 
a known voltage through the paper. 

The results obtained in this investigation compared 
favorably with those in Ref. 7. The ability to deter- 
mine if the boundary layer was turbulent or laminar became 
more difficult as the Reynolds number was increased and as 
the angle of incidence was increased. 
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II. DESCRIPTION OF APPARATUS 



A. LIQUID CRYSTALS 

Cholesteric liquid crystals are encapsulated in a poly 
vinyl alcohol binder which extends their useful life to 
several years. When applied to a black surface, the color 
of the liquid crystals will shift from red through the 
entire color spectrum to violet as the temperature is 
raised. The black surface ensures that all light passing 
through the liquid crystal film is absorbed, allowing the 
selectively scattered colors from the cholesteric liquid 
crystals to be observed. The process is reversible, 
repeatable, and can be accurately calibrated as to tempera 
ture within 0.1 degrees Celsius. 

Since each crystal has its own sensitivity and is 
colorless above and below its operating range, by select- 
ing the proper formulation of crystals a large number of 
isotherms can be viewed on a given surface at one time. 

An example would be the crystal S-40 (40 being the 
temperature in degrees Celsius at which the onset of red 
appears). As the temperature is raised to 40.0®C, the 
colorless crystals start to show red. The color will 
shift to green at 40.7°C, then to blue at 41.3°C, and 
finally to violet prior to becoming colorless. 
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Chloresteric liquid crystals are commercially available 
beginning at essentially any temperature from -20°C to 250°C. 
Reference 2 contains a very thorough background on the de- 
velopment of liquid crystals in its applications to the 
study of convective heat transfer. 

B . TEST BLADE 

The test blade profile used in this investigation is 
the NACA 65 Aj^q series blade with a span of ten inches and 
a chord of eight inches. It was originally constructed with 
an epoxy core covered by a thin metal casing. The blade is 
a typical airfoil shape used in compressors and was procured 
as original equipment for the Rectilinear Cascade at the 
U.S. Naval Postgraduate School Turbopropulsion Laboratory. 

The exact translation of the blade numbering system along 
with the two profile fairing and camber distribution equa- 
tions can be found in Ref. 8, 

C. LOW TURBULENCE SUBSONIC WIND TUNNEL 

The low turbulence subsonic wind tunnel is of an open 
circuit design with its intake inside the building and its 
exhaust outside. The test section is 20 inches by 28 
inches at the center and is 8 feet long. The prime mover 
is a six blade axial fan, located at the downstream end, 
driven by a variable speed 75 horsepower electric motor. 

The wind speed is continuously variable from 30 to 300 feet 
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per second when the test section is clear. When the blade 
and mounting apparatus is installed, the wind speed ranged 
from 10 to 290 feet per second. The freestream turbulence 
intensity is controlled by means of a six inch thick honey- 
comb flow straightener located at the plenum entrance, 
followed by up to five interchangeable graded screens and 
an area contraction of ten to one. 

Results of a complete survey of the test section are 
contained in Ref. 9, and indicated a turbulence in- 
tensity range from 0.05% to 0.20%. 
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III. EXPERIMENTAL PROCEDURE 



A. TEST BLADE PREPARATION 

Previous attempts to use this NACA 65 series airfoil 
for liquid crystal thermography were unsuccessful due to 
a short circuit between the blade and the covering material 
[Ref. 10] . This short circuit could not be eliminated with 
several coats of TV corona dope and red glyptol. To insure 
against such a problem during this investigation, the metal 
plating was completely removed. This plating measured 0.003 
inches thick. 

To better monitor the surface temperature, eight copper 
constantan thermocouples were imbedded in the blade along 
with the two previously installed copper constantan theinno- 
couples. The locations of the ten thermocouples are shown 
in Figure 1. The epoxy core was then covered with three 
coats of red glyptol. Prior to attaching to the blade, the 
temperature monitoring system consisting of the eight 
additional thermocouples, output meter, and connecting 
wires were calibrated using the Rosemo\int Engineering 
Company Model 913A calibration bath with distilled water 
and the Model 920A commutating bridge. The output meter 
used in this investigation is an Omega 2176A Digital 
Thermometer which reads in two-tenths of a degree Fahrenheit. 
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Figure 1- NACA 65 test blade showing location of thermo- 
couples, a) top view, b) side view. 
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The previous investigation [Ref. 10] has unsuccessfully 
attempted to use aluminized milar or Teledeltos paper as 
an electrically conducting cover material for the blade. 

A commercially produced electrically resistive paper known 
as Armstrong Temsheet was selected for this investigation 
due to its success reported in Ref. 5. Temsheet is 
a thin, highly flexible, electrically resistive, carbon 
impregnated paper, containing no wires or ribbons. The 
nominal thickness is 0.038 inches with an electrical re- 
sistance of approximately 45 ohms per square. The electri- 
cal resistivity is uniform, to within two percent, from 
point to point over large areas. Scotch Brand Sprament 
Adhesive was used to attach the temsheet to the blade. 

Self-adhering copper strips, attached to the edges of 
the temsheet, served as connecting points for the power 
supply. Power to heat the temsheet was provided by a 
Lambda model LK 34 5A FM regulated D.C. power supply. 

A Data Precision Model 1450 voltmeter, connected to the 
power supply, provided a digital readout of volts used to 
heat the temsheet. The resistance of the system was 
measured and found to be 11.7 ohms with no voltage and 
11.5 ohms with 40 volts. 

A mixture of encapsulated cholesteric liquid crystals, 
which consisted of R-33, S-40, R-45, and distilled water, 
was sprayed onto the center seven inches of the test blade 
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using an X-Acto magic touch air brush with 30 pounds per 
square inch air pressure. The distilled water was added 
in the amount necessary to allow the mixture to flow freely 
through the air brush. After spraying a single coating, 
the blade was subjected to a portable electric hot air 
blower to assist the drying of the temsheet and to check 
the color pattern of the crystals. With ten coats, the 
colors of the various crystals were adequate to provide the 
required temperature information. 

An investigation by Field in Ref . 5 noted that the 
liquid crystals fill-in irregularities in the surface of 
the temsheet leaving a roughness projection size of approxi- 
materly 0.00005 feet. 

The calibration results for the liquid crystals used 
are listed in Appendix B. Figure 2 shows the blade after 
the temsheet and copper strips were installed. 

B. BLADE INSTALLATION 

The NACA 65 series blade was attached to 8.5 inch dia- 
meter disks constructed from 0.5 inch thick plexiglass. 

A blade profile 0.25 inches deep was cut in each disk. 

The disks were then fitted to the blade and any gaps were 
filled in with a liquid epoxy to insure a vibration-free 
slip fit. The disks were fitted into the center of 20 inch 
by 29 inch plexiglass side panels, 0.5 inches thick, with 
a half lap joint to give a flush fit when they were mated. 
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Figure 2 . 



NACA 65 series test blade after temsheet and 
copper strip installation. 
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The blade was fitted with a handle so that rotating the 
blade and disks could be accomplished from outside of the 
wind tunnel. Rotating the blade could thus be accomplished 
without having to alter the wind speed controls to increase 
the angle of incidence. The disks were marked in 5 degree 
increments on the forward facing 180 degree arc, thereby 
allowing the angle of incidence to be easily observed 
through the clear side of the tunnel test section. The 
side panels were marked as to horizontal and vertical 
position. A plexiglass window was installed in the bottom 
of the test section, allowing viewing of the blade from 
underneath the tunnel. The test blade and side panels were 
rigidly mounted in the low turbulence wind tunnel above the 
viewing windown. The location of the isotherms on the 
blade were recorded with the aid of distance markings placed 
on each side of the temsheet and straight lines fixed to the 
viewing window. 

Figure 3 shows the final blade setup in the plexiglass 
support frame. Figure 4 shows the support frame, test 
blade, and electrical hookup inside the wind tunnel. 

C. BLADE SURFACE TEMPERATURE 

The surface temperature on the test blade was determined 
by the color of the liquid crystals. The thermocouples 
acted as an aid to determine which crystal was actually being 
viewed. Appendix B gives the temperatures corresponding to 
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Figure 3. Test blade and support frame prior to installa- 
tion in the low turbulence subsonic wind tunnel. 
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Figure 4. Support frame with the test blade and electrical 
hook-up inside the low turbulence subsonic wind 
tunnel . 
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the crystal/color combination. The location of the start 
of a certain color band was recorded as the line of the 
isotherm with the width of the color band giving an indi- 
cation of the temperature gradient along the blade. 

An increase or decrease of the applied voltage with 
the Reynolds number and angle of incidence held constant, 
raised or lowered, respectively, the overall surface 
temperature of the blade. This moved the color bands over 
the blade, making it possible to obtain a continuous tempera- 
ture profile of the entire blade by recording the location 
of the isotherm and crystal/color combination of these iso- 
therms. Scales were attached to each side of the blade, 
marked in 0.05 inch increments. Sighting lines on the 
viewing window enabled the isotherm to be recorded within 
a tolerance of 0.1 inch distance from the blade's leading 
edge . 

To insure a steady state condition when recording the 
blade surface temperatures, the thermocouples imbedded in 
the blade under the temsheet were monitored. Once the 
thermocouple nearest the desired isotherm did not fluctuate 
more than two-tenths of a degree in a three-minute period, the 
blade was considered to be in a steady state temperature condition. 

D. ANGLE OF INCIDENCE 

The actual tilt of the blade with respect to the free 
stream air direction was recorded as the angle of incidence. 
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This angle is defined, for this investigation, as the angle 
between the free stream air direction and a base line drawn 
between the lowest points on the leading and trailing edges 
of the blade. These baseline points were established by 
placing the test blade on a flat surface with the suction 
surface of the blade up, and marking the pressure surface 
of the blade at the points of contact between itself and 
the flat surface. Figure 5 shows a 15° angle of incidence. 
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IV. RESULTS 



The flow of a real fluid can be analyzed by the changes 
occurring within the boundary layer that is formed near the 
surface. Initially, the boundary layer developed is laminar. 
In laminar flow, fluid particles move along smooth paths 
with one layer gliding smoothly over an adjacent layer. 

At some critical distance from the leading edge, depending 
on the flow field and fluid properties, small disturbances 
in the flow begin to become amplified, and a transition 
process takes place until the flow becomes turbulent. The 
turbulent flow is characterized by a random churning action 
with chunks of fluid moving in irregular patterns within 
the boundary layer. 

The convection heat transfer coefficients and tempera- 
ture distributions are functions of both the thickness of 
the boundary layer and the type of flow within this 
boundary layer. As the boundary layer increases in thickness 
with distance from the leading edge and angle of incidence, 
the fluid velocity is retarded, resulting in an increased 
thennal resistance to heat transfer. This increased thermal 
resistance shows up in a decreased value of the convection 
heat transfer coefficient. The value of the convection 
heat transfer coefficient would then be smallest at a point 
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where the kinetic energy of the fluid is attenuated 
sufficiently so that the fluid in the boundary layer 
cannot overcomes the adverse pressure gradient present 
on the surface. At this point the boundary layer separates 
from the surface. 

Flow as classified by Roshko in Ref, 11, can be 
subcritical, supercritical, or transcritical . Subcritical 
flow contains only laminar separation. The supercritical 
flow is characterized by a laminar separation bubble, re- 
attachment, followed by turbulent separation. With the 
disappearance of a separation bubble, the flow is in the 
transcritical range and the separation is not purely 
turbulent. 

Figure 6 depicts a typical subcritical flow pattern 
and the resulting trends in the convection heat transfer 
coefficient for flow over a compressor blade. Figure 7 
depicts the supercritical flow pattern, with its resulting 
trends in the convection heat transfer coefficient for 
flow over a compressor blade. 

First, an example of a laminar boundary layer and 
separation is seen in Figure 8. The value of the local 
convection heat transfer coefficient continuously decreases 
from the leading edge to approximately 68% of the distance 
from the leading edge point. At this point the value of 
the coefficient increases at a relatively steeper slope. 
Knowing the angle of incidence (5 degrees) , and the Reynolds 
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Sep aration 




% DISTANCE SUCTION SIDE 



Figure 6. Schematic of a compressor blade in subcritical 

flow with the resulting convection heat transfer 
coefficient distribution. 
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Figure 7. Schematic of a compressor blade in supercritical 
flow with the resulting convection heat transfer 
coefficient distribution. 
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number (104,000 for this configuration), it is inferred 
that a laminar boundary layer exists until the 68% distance 
point is reached. The flow then separates from the blade 
with a turbulent wake remaining over the last 32% of the 
blade. Figure 9 is a color photograph showing the isotherms 
using the previously discussed liquid crystals. Starting 
at the leading edge, the crystals seen are R-33, S-40, and 
R~45, respectively, until the 68% distance point, then the 
reverse sequence is seen with much narrower color bands. 

Holding the Reynolds number essentially constant and 
increasing the angle of incidence to 10 degrees, a separa- 
tion bubble, as discussed in Ref* 11/ can be produced. 

With the test blade in this configuration, the bubble 
begins at the point of laminar separation, which now has 
moved closer to the leading edge. From Figure 10, this 
point of laminar separation and start of the separation 
bubble is located at the 35% distance point. Over the 
region of the bubble, the value of the heat transfer co- 
efficient increases due to the cooling effect of the re- 
circulating air inside the bubble [Ref. 12] . At the 70% 
distance point, the value of the heat transfer coefficient 
again starts to decrease. This indicates that the flow 
has now reattached to the blade and a turbulent boundary 
layer is present. At the 95% distance point, this 
boundary layer separates to form a turbulent wake as indi- 
cated by an upswing in the curve. Figure 11 shows the 
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Figure 9. NACA test blade, i = 5 degrees, Re^ = 104,000. 

Suction side, air flow from top to^bottom. 
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Figure 11. NACA test blade, i = 10 degrees. Re = 104,000. 

Suction side, air flow from top to ^bottom. 
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isotherms on the test blade at an angle of incidence of 
10 degrees and a Reynolds number of 104,000. The color 
bands have not yet come together at the 70% distance point 
which will indicate the cooler part of the blade and the 
point of flow reattachment. 

Again holding the Reynolds number relatively constant 
(Re^ = 104,000) and increasing the angle of incidence to 
25 degrees, a condition is observed that indicates that the 
blade has passed through its stall angle and the flow is 
completely separated from the suction side of the blade. 

As would be expected. Figure 12 shows the values of the 
heat transfer coefficient to be relatively low. The rise 
in their values in the mid section and trailing edge areas 
of the blade, indicate a cooling effect produced by the re- 
circulating flow in these areas. Figure 13 shows the 
pattern of isotherms which correspond to this situation. 

The observance of the sequences of colors for a particu- 
lar crystal will give the direction of the temperature 
gradient. This, along with the discrete temperature informa- 
tion available at any point over the blade surface, enables 
one to accurately deteimiine the location of the points of 
interest. Figure 14 compares the local convection heat 
transfer coefficient distribution for four different Reynolds 
numbers with the angle of incidence fixed at 15 degrees. 
Laminar separation occurred at 5% distance from the 
leading edge for the lowest Reynolds number of 104,000. 
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c 



Figure 13. NACA test blade, i = 25 degrees, Re = 104,000. 

Suction side, air flow from top to Bottom. 
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This point then moved forward to 2% distance for the 
Reynolds number of 208,000. For the Reynolds numbers of 
454,000 and 617,000 the location of the point of laminar 
separation, transition, and reattachment was the leading 
edge. Again for the two lower Reynolds numbers, the first 
peak in the curves show the reattachment positions. The 
turbulent boundary layer separation is at the point of 
lowest heat transfer coefficient value. As can be seen, 
this occurs from 75% to 85% distance depending on the 
Reynolds number. 

Finally, investigations were conducted at higher 
Reynolds numbers. By positioning the test blade at 0 
degree angle of incidence, increasing the tunnel wind 
velocity to 100 feet per second, and then rotating the 
blade to an angle of incidence of 25 degrees, gave a 
Reynolds number of 414,000. Data was taken in the normal 
manner between the 60% distance point and the trailing 
edge. This is shown in Figure 15. In the area from the 
leading edge back xintil the 60% point, no temperature data 
was possible due to the instability of the isotherms. The 
colors were very erratic, ranging from no color at all to 
a mixture of all colors in no noticable pattern. Figure 16 
is a picture of the test blade in this state, with no iso- 
therms visible in the forward 60% of the blade. It is 
inferred that the blade was then in a condition of stall 
in which the fluid was detaching and reattaching to the 
foirward 60% of the suction side of the blade. 
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Figure 16. NACA test blade, i = 25 degrees. Re = 414,000. 

Suction side, air flow from top to ^bottom. 
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The blade was then returned to an angle of 0 degree, 
and the air speed further increased to 150 feet per second. 
The blade was rotated to an angle of incidence of 5 degrees 
for a Reynolds number of 614,000. Figure 17 shows the iso- 
therms of the fluid with the blade in this condition. These 
isotherms can be seen at the leading edge and between the 
50% distance point to the trailing edge, and were of the 
usual pattern, perpendicular to the direction of flow. 

From just after the leading edge to the 50% distance point, 
isotherms were seen parallel to the direction of flow. 

This area should be observed with a flow visualization 
technique such as smoke to verify their actual existence 
on the suction side of the test blade. The actual presence 
of such a complex three-dimensional flow could be the 
subject of further studies. 
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Figure 17. NACA test blade, i = 5 degrees, Re = 614,000. 

Suction side, air flow from top to bottom. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The liquid crystal thermography studied in this investi- 
gation shows that it's use will provide an excellent technique 
to obtain both quantitative and qualitative heat transfer 
information on heated objects placed in a forced convective 
environment. It was possible to quickly reduce the data to 
obtain information as to the convection heat transfer co- 
efficient as it varied over the entire suction side of the 
compressor blade. From the values of the heat transfer co- 
efficient, the character of the flow was easily determined. 
The location of separation point, reattachment point, 
separation bubbles, and stall angle were able to be deter- 
mined. With adjustments in the voltage supplied to the 
temsheet, a continuous picture was seen, which is very 
advantageous as compared to pressure taps or theirmocouples , 
which provide information only at discrete points. This 
continuous coverage enabled the viewing of the isotherms 
that were produced parallel to the direction of flow, 
particularly on the pressure side of the blade. 

The ability to observe the isotherms could be made 
easier if blades in the future would be mounted in a 
vertical position, with a viewing window to see both sides 
of the blade. To study the flow over a particular blade 
shape, the thickness of the temsheet would have to be 
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considered. A blade built to allow the temsheet to be 
installed in a recessed area with the surface of the temsheet 
conforming to the blade specifications, made from a material 
with a known value of the thermal conductivity, would improve 
the usefulness of the data. An accurate value of the thermal 
conductivity of temsheet should also be deteirmined, if this 
material is to be used for future investigations. If informa 
tion is needed at higher Reynolds numbers or at a greater 
angle of incidence, a low numbered crystal should be used, 
due to the voltage limitations of the equipment. 

As a final recommendation, prior to investigating the 
flow over the compressor blade in a cascade wind tunnel, 
the type of flow inferred from the convection heat transfer 
coefficient should be viewed with some type of flow visualize 
tion technique in the subsonic wind tunnel. 
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APPENDIX A 



ELECTRICAL POWER REQUIREMENT FOR LIQUID CRYSTAL THERMOGRAPHY 



Turbulent flow of air is assumed over a "flat plate" 
with a span of ten inches and a chord length of eight 
inches. Surface area equals 160 square inches. Plate 
surface temperature to be 100 degrees Fahrenheit. Air 
flow equals 275 feet per second at 65 degrees Fahrenheit. 

T ^ ^b+^ ^ l.Q.Q+$5. = 82.5°F 

f 2 2 



Air properties at 82.5°F: 

Ibm 



y = 0.04797 



hr-ft 



k 



0.0152 



Btu 

hr-ft-'^F 



Pr = 0.708 



p = 0.0735 



Ibm 
ft 3 



V = 1.8126 X 10 



-4 ft^ 
sec 






The above conditions give: 



Re = = 1.011 X 10^ 

c y 
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For turbulent flow, to find the average convective heat 



transfer coefficient, equation 5-81 of Ref* 13 is used: 



Nu^ = 



= Pr‘/> 
k 



0 • a 



Pr^/ " (0.037Re^ - 850) 

i-j 



hL 



= 1341.85 



h = 30.59 



Btu 

hr-ft2-°F 



Equating the heat generated in the test blade by the 
electrical current to the heat removed by forced convection 
(heat loss by radiation was considered negligible) gives: 



Power = I^R = hA(Tj^ - T„) = VI 



which leads to 



I . 
min 



hA(Tj^ - T^) 
V 

max 



P 



V 

max 



and 



R 



max 



hA(Tj^-T^) 



12 



min 



P 



12 . 
min 



The power required to heat the blade is calculated as 
P = hA(T, - T ) = 349 watts 

D °° 

For a power supply rated at a maximum output potential 
of 60 volts, current and resistance requirements are : 

I . =5.8 amps 

min 

R =10.3 ohms 
max 
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APPENDIX B 



LIQUID CRYSTAL CALIBRATION RESULTS 



Liquid Crystal 




Color Transition 




(NCR Designation) 


Red (°F) 


Green (°F) 


Blue i° 


R-33 


91.4 


92.5 


94.2 


S-40 


104.1 


105.2 


106.3 


R-45 


106.6 


108.2 


109.5 
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APPENDIX C 



BLADE EDGE EFFECTS 



A. SIDE WALLS 

To eliminate any readings that could be affected by the 
boundary layer caused by flow pass the two side walls, the 
boundary layer thickness was calculated for a "worst case" 
situation. 

Assume fully turbulent boundary layer, equation 5-85 of 
Ref. 13 is used; 



6 

X 



0.381 Re~^/^ 

X 



which leads to 



6 = 0.381X''/® (Puoo)-i/5 

y 






6 = 0.03ix''^^ inches 



for X (distance from leading edge of side wall to 
trailing edge of blade) equal to 16 inches 



6 = 0.286 inches 



Re 



X 



368,000 
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B. TEST BLADE 



The span wise variation in the temperature field due 
to end losses was estimated by treating the temsheet as a 
thin fin. Symmetry about the centerline was assumed. A 
"worst case" situation was also assumed with the temperature 
at the outer edge of the temsheet equal to the free stream 
temperature, T^. Using the equation developed in Ref. 5: 

T-Tqo _ Cosh (my) . 

T -T ^ Cosh(mL)^ 
o '=° 

where 

y = distance from centerline 
L = distance from centerline to outer edge 
T = temperature at location y 
(T^ - Too) = centerline temperature excess 




t = thickness of test section material 
k = thermal conductivity of test section material 
h = surface heat transfer coefficient 
T^ = free stream air temperature 
The boundary conditions used for the problem were 




at y = L, T = T, 
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Using representative values of: 
Btu 



h = 5 



hr-ft^-®? 



k = 0.06 



Btu 



hr-ft-°F 
t = 0.039 inches 



[Ref. 5] 



[Ref. 5] 



2L = 8 inches 

The distance y for which ( 



T-Tg. , _ 



To-T. 



) = 0.99 was found to 



be 3.65 inches. Figure 18 shows the final placement of the 
temsheet and copper strips on the test blade. 

It was concluded that if the electrodes were spaces 
8 inches apart and 1 inch from the ends of the blade, using 
only the center 7 inches, would yield temperature informa- 
tion that was essentially uninfluenced by effects from the 
turbulent boundary layer developed by the side walls and 
heat lost through the edges of the temsheet. 
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i span ► 







Figure 18. 



Placement of temsheet and copper strips 
(electrodes) on the test blade. 
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APPENDIX D 



DATA AND DATA REDUCTION 



In order to experimentally determine the convection 
heat transfer coefficient as a function of location on the 
suction side of the blade, measured from the leading edge, 
a relationship between the surface heat flux produced by 
the Joulean heating effect in the temsheet and the losses 
through radiation and conduction had to be established. 
Performing a simple heat balance on the blade gives 



h 



^electrical ^conduction 



h 



c 




r 



where 



'^electrical 



= Heat transfer into blade 




= Heat transfer loss through the two ends 



not covered by the temsheet 



h = Radiation heat transfer coefficient 
r 



= [oF^_2(T+TJ (T^+T^)] 

F, „ = Radiation exchange factor between the test blade 
and its surroundings 



a = Stefan-Boltzman constant 



A 



surface area of the short end fins 



s 
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The electrical heat transfer value was determined by 



Q 



electrical 



Yl 

R 



The conduction losses through the two blade ends were 

estimated by modeling the ends as short fins, with convection 

end losses. A value of h was estimated using the average of 

the h calculated without the conduction loss factor for the 
c 

Reynolds number and angle of incidence of the blade. Using 
case "2" from Ref. 13 (fin of finite length and losses) 
by convection from its end) 



(JL) 

n = fh^R ^kA Mmx sinh(mL) + (mk) cosh(mL) , 

^conduction P cosh(mL) + (J^) sinh(mL) s 

mk 

where 

h = average convection heat transfer coefficient for the 
short end fins 

B = blade profile perimeter 
P 

k = thermal conductivity of the test blade 

A = profile area of the test blade 
P 

AT = T-T„ 



hB 






L = length of the short fin 

A = surface area of the short fin 
s 
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A sample calculation is provided for illustration using 
the following values from the data at a Reynolds number of 



208,000, a 15 degree angle of incidence, and 7.0 inches 
from the leading edge on the suction side of the blade: 
V = 31.45 volts 
R = 11.5 ohms 



B_ = 16.1 inches 
P 



k = 0.2 



Btu 

hr-ft-^’F 



T = 109. 5°F 

Too = 60.8°F 

A =4.6 in^ 
P 



A^ = 16.6 in' 
s 



L = 1 inch 
t = 0.75 inch 
a = 0.1714 X 10 



-8 Btu 



hr-ft2-°R^ 



Fi _2 = 1.0 



h = 16 



Btu 



hr-ft2-°F 
mL = 16.99 
h = 1.11 



Btu 



hr-ft2-°F 



^electrical = ”3.46 ^ 



Q . 2.11 (for each end) 

conduction = hr 



Btu 



Btu 



hr-ft2-°F ■ 0-‘752 hr-ft2-°F 



- 1.11 



Btu 



hr-ft2°F 



. _ 4 71 Btu 

hr-ft2=°F 
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The data sheets for all four Reynolds numbers investi- 
gated at an angle of incidence of 15 degrees follow. 
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Air Speed 0.25 in H 2 O 
25 ft/sec 



DATA 



Re_ 104,000 



i angle 15 degrees 



h 



10.0 



Btu 

hr-ft2"^^ 



P ^ 30.07 in Hg 

atm 



V 


R 


Cry. 


Color 


Tenp. 

(°F) 


Location 
(in) (%) 


Too 


AT 


Uncorr 

h 

c 


h 

r 


Corr 

h 

c 


27.77 


11.6 


R-33 


R 


91.4 


5.3 


63.5 


63.0 


28.4 


8.71 


1.06 


7.05 


27.77 


11.6 


R-33 


G 


92.5 


5.5 


65.9 


63.0 


29.5 


8.39 


1.07 


6.73 


27.77 


11.6 


R-33 


B 


94.2 


5.6 


67.1 


63.0 


31.2 


7.93 


1.07 


6.26 


27.77 


11.6 


S-40 


R 


104.1 


6.0 


71.9 


63.0 


41.1 


6.02 


1.10 


4.32 


27.77 


11.6 


S-40 


G 


105.2 


6.1 


73.1 


63.0 


42.2 


5.86 


1.11 


4.16 


27.77 


11.6 


S-40 


B 


106.3 


6.2 


74.3 


63.0 


43.3 


5.71 


1.11 


4.01 


27.77 


11.6 


R-45 


R 


106.6 


6.5 


77.8 


63.0 


43.6 


5.68 


1.11 


4.04 


27.77 


11.6 


R-45 


G 


108.2 


6.75 


80.8 


63.0 


45.2 


5.47 


1.12 


3.76 


27.77 


11.6 


R-45 


B 


109.5 


7.25 


86.8 


63.0 


46.5 


5.32 


1.12 


3.61 


27.77 


11.6 


R-45 


G/R 


107.4 


7.75 


92.8 


63.0 


44.4 


5.57 


1.11 


3.87 


27.77 


11.6 


S-40 


B/G/R 


105.2 


8.0 


95.8 


63.0 


42.2 


5.86 


1.11 


4.16 


27.77 


11.6 


R-33 


B/G/R 


92.8 


8.15 


97.6 


63.0 


29.8 


8.30 


1.07 


6.64 


36.85 


11.6 


R-33 


R 


91.4 


1.25 


15.0 


63.0 


28.4 


15.34 


1.06 


13.68 


36.85 


11.6 


R-33 


G 


92.5 


1.35 


16.2 


63.0 


29.5 


14.77 


1.07 


13.11 


36.85 


11.6 


R-33 


B 


94.2 


1.5 


18.0 


63.0 


31.2 


13.97 


1.07 


12.30 


36.85 


11.6 


S-40 


C 


105.2 


4.5 


53.9 


63.0 


42.2 


10.32 


1.11 


8.62 


36.85 


11.6 


R-33 


B/G/R 


92.8 


0.5 


6.0 


63.0 


29.8 


14.62 


1.07 


12.96 


36.85 


11.6 


R-33 


R/G/B 


92.8 


0.2 


2.4 


63.0 


29.8 


14.62 


1.07 


12.96 
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Air Speed 0 . 5 in H 2 O 
50 ft/sec 

Re^ 208,000 



DATA 



1 Angle 15 degrees 

Btu 



h 15.0 



hr-ft2-°F 



24.99 in Hg 



V 


R 


Cry. 


color 


Tertp. 

(°F) 


Location 
(in) (%) 


Too 


AT 


Uncorr 


hr 


Corr 

he 


31.45 


11.5 


R-33 


R 


91.4 


5.4 


64.7 


60.8 


30.6 


10.46 


1.06 


8.68 


31.45 


11.5 


R-33 


G 


92.5 


5.5 


65.9 


60.8 


31.7 


10.10 


1.06 


8.31 


31.45 


11.5 


R-33 


B 


94.2 


5.6 


67.1 


60.8 


33.4 


9.58 


1.07 


7.79 


31.45 


11.5 


S-40 


R/G 


104.7 


5.8 


69.5 


60.8 


43.9 


7.29 


1.10 


5.47 


31.45 


U.5 


S-40 


G/B 


105.8 


5.9 


70.7 


60.8 


45.0 


7.11 


1.10 


5.28 


31.45 


11.5 


R-45 


R 


106.6 


6.2 


74.3 


60.8 


45.8 


6.99 


1.10 


5.16 


31.45 


11.5 


R-45 


G 


108.2 


6.45 


77.2 


60.8 


47.4 


6.75 


1.11 


4.92 


31.45 


11.5 


R-45 


B 


109.5 


7.0 


83.8 


60.8 


48.7 


6.57 


1.11 


4.73 


31.45 


11.5 


R-45 


G 


108.2 


7.25 


86.8 


60.8 


47.4 


6.75 


1.11 


4.92 


31.45 


11.5 


R-45 


R 


106.6 


7.50 


89.8 


60.8 


45.8 


6.99 


1.10 


5.16 


31.45 


11.5 


S-40 


B/G/R 


105.2 


7.75 


92.8 


60.8 


44.4 


7.21 


1.10 


5.38 


31.45 


11.5 


R-33 


B/G/R 


92.8 


8.0 


95.8 


60.8 


32.0 


10.0 


1.06 


8.21 


44.31 


11.5 


R-33 


R 


91.4 


1.15 


13.8 


60.6 


30.8 


20.63 


1.06 


18.85 


44.31 


11.5 


R-33 


G 


92.5 


1.3 


15.6 


60.6 


31.9 


19.92 


1.06 


18.13 


44.31 


11.5 


R-33 


B 


94.2 


1.6 


19.2 


60.8 


33.6 


18.91 


1.07 


17.12 


51.7 


11.5 


R-33 


R/G 


92.0 


0.5 


6.0 


60.8 


31.2 


27.73 


1.06 


25.94 


51.7 


11.5 


R-33 


G/B 


93.4 


0.65 


7.8 


60.8 


32.6 


26.54 


1.06 


24.75 


51.7 


11.5 


R-33 


R/G/B 


92.8 


0.05 


0.1 


60.8 


32.0 


27.03 


1.06 


25.24 


51.7 


11.5 


R-33 


B/G/R 


92.8 


0.25 


3.0 


60.8 


32.0 


27.03 


1.06 


25.24 
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Air Speed 2 . 6 in H2O 
110 ft/sec 



DATA 



Re^ 454,000 



i Angle 15 degrees 
Btu 



h 25 
^atrn 



hr-ft^-°F 
30.09 in Hg 



V 


R 


Cry. 


Color 


Tenp. 

(°F) 


Location 
(in) {%) 


Too 


AT 


Uncorr 


hr 


Corr 

he 


36.07 


11.7 


R-33 


R 


91.4 


4.75 


56.9 


67.6 


23.8 


17.39 


1.08 


15.37 


36.07 


11.7 


R-33 


G 


92.5 


5.0 


59.9 


67.6 


24.9 


16.62 


1.08 


14.60 


36.07 


11.7 


R-33 


B 


94.2 


5.2 


62.3 


67.6 


26.6 


15.56 


1.09 


13.53 


36.07 


11.7 


S-40 


P/G 


104.7 


5.5 


65.9 


67.6 


37.1 


11.16 


1.12 


9.10 


36.07 


11.7 


S-40 


G/B 


105.8 


5.6 


67.1 


67.6 


38.2 


10.83 


1.12 


8.77 


36.07 


11.7 


R-45 


R 


106.6 


5.75 


68.9 


67.6 


39.0 


10.61 


1.12 


8.55 


36.07 


11.7 


R-45 


G 


108.2 


5.95 


71.3 


67.6 


40.6 


10.19 


1.13 


8.13 


36.07 


11.7 


R-45 


B 


109.5 


6.3 


75.4 


67.6 


41.9 


9.88 


1.13 


7.80 


36.07 


11.7 


R-45 


G 


108.2 


6.65 


79.6 


67.6 


40.6 


10.19 


1.13 


8.13 


36.07 


11.7 


R-45 


R 


106.6 


6.9 


82.6 


67.6 


39.0 


10.61 


1.12 


8.55 


36.07 


11.7 


S-40 


B/G 


105.8 


7.25 


86.8 


67.6 


38.2 


10.83 


1.12 


8.77 


36.07 


11.7 


S-40 


G/R 


104.7 


7.35 


88.0 


67.6 


37.1 


11.16 


1.12 


9.10 


36.07 


11.7 


R-33 


B/G 


93.4 


7.8 


93.4 


67.6 


25.8 


16.04 


1.08 


14.02 


36.07 


11.7 


R-33 


G/R 


92.0 


8.0 


95.8 


67.6 


24.4 


16.96 


1.09 


12.65 


56.7 


11.7 


R-33 


R 


91.7 


0.45 


5.4 


68.2 


23.8 


44.08 


1.08 


42.06 


56.7 


11.7 


R-33 


G 


92.5 


0.65 


7.8 


68.2 


24.3 


42.09 


1.08 


40.07 


56.7 


11.7 


R-33 


B 


94.2 


1.0 


12.0 


68.2 


26.0 


39.34 


1.09 


37.31 


56.7 


11.7 


R-45 


G 


108.2 


3.9 


46.7 


68.2 


40.0 


25.57 


1.13 


23.50 
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DATA 



Air 


Speed 5.0 


in H^O 


i Angle 15 degrees 




150 


ft/sec 


Btu 

^ — hr-ft2-°F 


Re 

c ■ 


617,000 




P ^ 30.02 in Hg 

atm ^ 



V 


R 


Cry. 


Color 


Toip. 

(°F) 


Location 
(in) (%) 


Too 


AT 


tncx)rr 

h 

c 




Corr 

he 


43.1 


11.5 


R-33 


R 


91.4 


4.5 


53.9 


67.6 


23.8 


25.26 


1.08 


23.15 


43.1 


11.5 


R-33 


G 


92.5 


4.65 


55.7 


67.6 


24.9 


24.15 


1.08 


22.04 


43.1 


11.5 


R-33 


B 


94.2 


4.8 


57.5 


67.6 


26.6 


22.60 


1.09 


20.49 


43.1 


11.5 


S-40 


R/G 


104.7 


5.25 


62.9 


67.6 


37.1 


16.21 


1.12 


14.06 


43.1 


11.5 


S-40 


G/R 


105.8 


5.3 


63.5 


67.6 


38.2 


15.74 


1.12 


13.59 


43.1 


11.5 


R-45 


R 


106.6 


5.6 


67.1 


67.6 


39.0 


15.42 


1.12 


13.26 


43.1 


11.5 


R-45 


G 


108.2 


5.7 


68.3 


67.6 


40.6 


14.81 


1.13 


12.65 


43.1 


11.5 


R-45 


B 


109.5 


6.35 


76.0 


67.6 


41.9 


14.35 


1.13 


12.19 


43.1 


11.5 


R-45 


G 


108.2 


6.75 


80.8 


67.6 


40.1 


14.81 


1.13 


12.65 


43.1 


11.5 


R-45 


R 


106.6 


7.0 


83.8 


67.6 


39.0 


15.42 


1.12 


13.26 


43.1 


11.5 


S-40 


B/G/R 


105.2 


7.5 


89.8 


67.6 


37.6 


15.99 


1.12 


13.84 


43.1 


11.5 


R-33 


B/G/R 


92.8 


8.15 


97.6 


67.6 


25.2 


23.86 


1.08 


21.75 


57.71 


11.5 


R-33 


G 


92.5 


1.5 


18.0 


68.2 


24.3 


44.36 


1.08 


42.25 


57.71 


11.5 


R-33 


B 


94.2 


1.85 


22.2 


68.2 


26.0 


41.45 


1.09 


39.34 
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APPENDIX E 



UNCERTAINTY ANALYSIS 



The uncertainties for the final results were calculated 
using the method of Kline and McClintock described in 
Ref. 14. The measured variables which were the origins 
of uncertainty were voltage, resistance, free stream air 
temperature, blade surface temperature, free stream air 
velocity, atmospheric conditions, angle of incidence, 
isotherm location, and blade physical measurements. 



THE REYNOLDS NUMBER 

The Reynolds numbers were calculated from 



Re = 
c 



uc 

V 



A Re 




V 



+ 




+ 



(^) 

u 




For the experimental run at a free stream velocity of 
25 feet per second 

.-4 ft^ 



V = 1.6 X 10 + 0.1 X 10 

c = 0.67 + 0.01 ft 
ft 



sec 



u = 25 + 10% 



sec 
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ARe 
c 

Rec 



^ 0.1 X 10 
. 1.6 X 10 



•) + ( 



0.01 

0.67 



2 

) + ( 



.5 

25 




ARe 

c 



Re 



c 



0.119 



with Re = 105,000 
c 



ARe^ = 12,000 



The values of the uncertainties for the results are 
listed below for a Reynolds number of 105,000 and an angle 
of incidents of 0 degrees. 



Quantity 

hr 

^conduction 

*^electrical 



Uncertainty 

9.2% 

17.7% 



0.9% 

20 . 0 % 



The values of the uncertainties for the measured 
variables are listed below. 



Quantity 


Uncertainty 


V 


0.01 volts 


R 


0 . 1 ohms 


OO 

X (isotherm 
location) 


0.2°F 
0 . 1 inch 


Tfc X 


0.1°F 


i 


1 . 0 degree 


^atm 


0.01 in. Hg 


u 


10% 
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